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Quantum breathing mode of trapped bosons and fermions at arbitrary coupling 
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Interacting particles in a harmonic trap are known to possess a radial collective oscillation - the 
breathing mode (BM). We show that a quantum system has two BMs and analyze their properties 
by exactly solving the time-dependent Schrodinger equation. We report that the frequency of one 
BM changes with system dimensionality, the particle spin and the strength of the pair interaction 
and propose a scheme that gives direct access to key properties of trapped particles, including their 
many-body effects. 
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I. INTRODUCTION 

The dynamics of trapped quantum systems are of 
growing interest in many fields, including c orrela ted elec- 
trons in metal cluster»2 or quantum dot^^IUll and ul- 
tracold Bose and Fermi gases in traps or optical lat- 
tices, for recent overviews see e.g.l^. Particular atten- 
tion has recently been devoted to Bose-Einstein conden- 
sation in low dimension and to the analysis of nonide- 
ality (interaction) eSect^^^^^^^, including superfluidity 
and crystallization, se^^. At the same time, nonideal 
low-dimensional Bose and Fermi systems present major 
experimental and theoretical challenges. For a reliable 
diagnostics of static and time-dependent propertie s, the 
collective oscillations of the system play a key rol d^^^^^l 
It is the purpose of this paper to show that, among them, 
the monopole or breathing mode (BM) which is easily ex- 
cited experimentallj0 carries particularly valuable infor- 
mation on the system dimensionality d, the spin statistics 
of the particles and on the form and relative strength A 
of their pair interaction given by Eqs. (2|3l. 

The BM of a system of particles in a harmonic trap 
with frequency is well known in two limiting cases: in 
the limit of very strong interparticle repulsion (A — *■ oo) 
with an inverse power law potential, wdrj) ^ r^\ the 
particles behave classically being spatially localized, and 
the BM which describes the radial expansion/contraction 
of a c loud o f well-separated particles is independent of N 
and (P^^, whereas its frequency is sensitive to the in- 
teraction, wbm = \/2r+Tf2. In the second limiting case, 
that of an ideal quantum gas with A = 0, the BM corre- 
sponds to a periodic expansion/contraction of the wave 
function with frequency wbm = 2f2, which is, of course, 
independent of N, d and the type of the interaction, see 
SecEEl 

This raises the question about lubm for arbitrary finite 
values of the interaction strength which we answer in this 
paper. While our results for the BM are representative 
for any nonideal quantum system in a harmonic trap, 
here, we concentrate on a complete analysis of the case of 
Coulomb interaction (I — 1). We present exact numerical 
results from solutions of the time-dependent Schrodinger 
equation (TDSE) for two fermions and bosons in Id 



and 2cf^ and time-dependent Hartree-Fock (TDHF) re- 
sults for A = 2 — 4. We show that, in fact, there co- 
exist two independent breathing modes one of which is 
A— dependent and the relative spectral weight of which 
varies with A. Further, fermions and bosons have the 
same breathing modes in Id, whereas a substantial dif- 
ference is observed in higher dimensions. 



II. MODEL AND SOLUTION PROCEDURE 

We start by considering two identical particles (mi = 
7712 = "T-i 91 — <l2 = q) in a harmonic potential of 
frequency with Coulomb repulsion, described by the 
TDSE of two particles 
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where the total potential V is the sum of harmonic con- 
finement and Coulomb repulsion. The coupling parame- 
ter ~ the ratio of the mean interaction and single-particle 
energy - is given by 



A = 



1 



(3) 



with the oscillator length Iq = ^hjmO,. Throughout this 
work, lengths, times and energies will be given in units 
of Zoi^^^^ and h£l, respectively. 



A. Initial conditions 

Eq. ([T]) has to be supplemented by an initial condi- 
tion for the two-particle wave function. Here we use ei- 
ther a symmetric (JS) or anti-symmetric (A) function, 
^S.A ^ 0) = ±^''^"^(r2,ri,0), the symmetry 

of which is preserved during the time evolution since the 
hamiltonian is spin-independent. For the actual form of 
^'q'^ we choose a stationary solution of Eq. ([ij which 
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FIG. 1: (color online) Time evolution of the one-particle 
potential energy Af/pot = (t/pot)(t) — Upot,t=o and A|x| = 
{\x\){t) — |a;t=o| for N = 2 particles in a ID trap at A = 1.0, 
obtained by solving Eq. |l|, symbols. Lines: fit, Eq. (j4|, with 
frequencies cur = 1.901 and ljr — 2.0. 

is computed by imaginary time stepping. To single out 
time-dependent solutions of Eq. (jlj of pure BM-type we 
use two different excitation methods: (I) a fast switch 
of the trap frequency and (II), the response to a pe- 
riodic modulation of fl^^. While the former has com- 
putational advantages, the latter is more easily realized 
experimentallj^^. 

B. Numerical solution of the TDSE 

We solve the two-particle TDSE ([ij numerically us- 
ing two independent methods: i) a standard grid-based 
Crank-Nicolson scheme (cn) with at least 1000 grid 
points in each direction and ii) by expanding the wave 
function into a basis of oscillator eigenfunctions (OB) 
using up to Ni, = 625 basis functions. The numerical 
parameters are adjusted for each value of A, such that 
the results are fully converged with respect to the time 
step size, the simulation box size and the number of grid 
points/basis functions. This has to be undertaken very 
carefully especially in the case of large A, where the re- 
pulsive Coulomb interaction attains large values and the 
accurate numerical treatment becomes challenging. 

III. RESULTS 

A. Dynamics following a switch of the confinement 

Let us start with an antisymmetric initial state "^/q and 
apply method (I) . This is realized by turning off the trap 
for a short time, typically At = 0.1 after which it is 
restored (the explicit time-dependence is not relevant as 
long as the excitation is spectrally much broader than il). 
During the "off" cycle, the Coulomb repulsion drives the 
particles out of their initial equilibrium state and initiates 
the BM. 

To quantify the oscillatory motion we compute the 
time-dependent expectation value of the single-particle 
potential energy (C/pot)(^) with respect to the external 



trap and, independently, the expectation value (|a:;|)(i) 
which is directly associated with the monopole oscilla- 
tion. Our numerical results confirm that both quantities 
exhibit identical time dependencies. In Fig. [TJwe show 
the result for = 2 in a ID trap at an intermediate 
coupling, A = 1, where we expect a quantum BM (QBM) 
with frequency inbetween the ideal quantum and classical 
limits, 2 and -v/S, respectively. However, the simulations 
reveal a different behavior with evidence of a beating of 
two oscillations. This is readily confirmed by applying a 
two-frequency fit 

f{u!r,ujR,t) ^ a-sm[uJr{t-to)]+b-sm[ujR{t-tQ)]+fo, (4) 

to our data (/ stands for (C/pot) or (|a;|)), where to and 
t'o indicate phase shifts, a and b are the amplitudes and 
ujr and lor the frequencies of the two QBMs, and /o is 
the unperturbed (equilibrium) value of /. With the val- 
ues ujR = 2 and LUr — 1.901 perfect agreement with the 
simulations is achieved, cf. lines in Fig. [l] 

The origin of the two QBMs becomes immediately 
clear from the structure of Eq. ([l]) . It can be solved with 
a product ansatz, ^'(R, r,i) = 0(R, i) •(/^(r, i), factorizing 
into functions of the center of mass (CoM) and relative 
coordinates, R = (ri + r2)/2, and r = ri — r2, resulting 
in two independent TDSEs: 

^dMr,t)=(-^ + \r' + ^y{r,t), (5) 

whereas 0(R, t) obeys a simple harmonic oscillator prob- 
lem which is independent of A: 

idtm,t)= (^-\^+-R.^^m,t) ■ (6) 

This well-known splitting indicates the existence of two 
independent motions related to the relative and the CoM 
problem corresponding, in case of excitation (I), to two 
QBMs with the characteristic frequencies LOr and ur, re- 
spectively. Obviously, Eq. ([6| leads, in our case, to the 
(ideal) QBM with frequency ujr = 2 because it is just the 
TDSE for a single harmonic oscillator with mass m = 2, 
se^lSI. Thus, for all couplings A, the system ^ possesses 
two QBMs, one with a universal frequency ujr — 2 and 
one with a A— dependent frequency LUr', the two modes 
are depicted by the arrows in the left part of Fig. |2] for 
an illustration, see Ref.^^. In the ideal quantum limit, 
A ^ 0, Eq. (|5| transforms to Eq. ([6| leading to a two- 
fold degenerate mode with frequency ioj^^. 

B. Resonance excitation. Absorption spectrum 

We now verify these results by applying the reso- 
nance excitation (II), i.e. using a periodically modu- 
lated trapping potential of the form ri^[l -I- P{t)]rf/2, 
with p{t) = (3q exp[— (t — t^Y /2a\ sin(a;oxti) and i = 1,2. 
Since we expect a strong resonance for Woxt = and 
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FIG. 2: Resonance excitation (II) for A = 1.0 and three fre- 
quencies, LJext = a(^r with Q = 0.95,0.995,1.0. Left: anti- 
symmetric initial state (top) and particle density (grey area). 
Arrows indicate the two breathing motions. Right: time evo- 
lution of the total energy for the three external frequencies 
(top) and the exciting pulse (bottom). 



cjcxt = we use a small modulation depth and a fi- 
nite pulse width, typically Pq — 5 ■ 10~^ , to = 240 and 
a = 100. The system response is characterized by the 
behavior of the total energy, E^tot, see Fig. |2]for typical 
examples. For a quantitative analysis we use the value 
E^^. sufficiently long after the excitation (cf. the dotted 
horizontal lines) . 

We now construct the absorption spectra by record- 
ing the values E^^{LOey,t, X), by scanning the frequency 
Wext for a given A and repeating this procedure for a 
broad range of A values. The resonance peaks are shown 
in Fig. |3] and confirm the existence of two frequencies. 
Their values fully agree with the result of method (I) , cf . 
the grey line. Moreover, from the peak areas we can de- 
duce the relative spectral weight of the two modes, clearly 
showing the continuous transition from the ideal quan- 
tum limit (where the mode ljr dominates) to the strongly 
coupled classical case. Both peaks merge at weak cou- 
pling, around A — 0.01 (for the given a). On the other 
hand, even for A as large as 20, the pure quantum mode 
is clearly observable, being of the same absolute im- 
portance as for A = 0. 

Furthermore, we underline that our approach of solv- 
ing the full two-particle TDSE is applicable beyond the 
linear response regime. It allows us to estimate the ex- 
act energy absorption spectrum including the line shapes 
which obviously depend on the type of the external per- 
tubation used. This is of special importance for experi- 
mental situations, where every excitation has finite du- 
ration and typical modulations which will also result in 
characteristic line shapes. 



C. Analytical approximations for cUr 

In addition to the numerical results it is possible to 
derive semi-analytical expressions for the breathing fre- 
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FIG. 3: QBM frequencies cjr and uir for N = 2 particles 
in a ID trap versus A from solution of Eq. ([5| with a basis 
expansion method (grey line) and a cn-solution of Eq. ([l]), (-I-) 
using excitation method (I). Grey resonance spectra (-EtS) are 
the result of method (II). Dashed lines: classical (A > 30) and 
quantum (A < 1) mean-field models, symbols: TDHF results 
for A'' = 2, 3, 4, which coincide for A < 1. 

quency LUr in the two limits of weak and strong coupling, 
respectively. To obtain uj^ for small A, we use a quan- 
tum mechanical Hartree model. Since for A 1, the 
interaction potential w is a small perturbation it can be 
approximated by the ansatz 

wir,~r,)^^-J2 ld\w{r,^r)\q>t'ir)\', (7) 

involving the densities of the ideal (undisturbed) one- 
particle states. This potential, together with the external 
potential, gives rise to a renormalized confinement with 
the effective harmonic trap frequency fi^*^ yielding the 
QBM frequency uir = 2ri°*^. As can be seen in in Fig. [s] 
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equation for iV > 3 with a sufficient accuracy is computa- 
tionally very expensive. Therefore, we performed instead 
time-dependent Hartree-Fock calculations for three and 
four particles using a harmonic oscillator representation 
of the HF equations^^, see the symbols in Fig. [s] Full 
agreement with the TDSE results is observed for N — 2, 
up to A ~ 3. Further, the data for = 3, 4 are very close 
to the TDSE results for = 2 up to A 1. From this we 
expect that also the exact Ci.'r(A) is only very weakly A^- 
dependent as in the ideal quantum (A = 0) and classical 
limitd^. However, a conclusive answer to this question 
requires a thorough computational analysis on the basis 
of the TDSE for larger particle numbers. 



FIG. 4: (color online) Frequency cUr for 2 particles with a 
symmetric initial state in a ID trap from numerical solution 
of the TDSE using a softened Coulomb potential. Crosses 
correspond to solution method i) and lines to method ii), see 
Section fll B| For k — 0, uir is the same for symmetric and anti- 
symmetric states. The result for the anti-symmetric case is 
shown by the monotonic full (black) line. A finite k drastically 
influences Wr for the symmetric case which is explained by the 
behavior of the relative wave function at a; = (see inset). 



the result is very close to the exact numerical solution of 
the TDSE for A < 0.3, cf. the dotted line. 

Similarly, for large A, we construct a semi-classical 
mean-field theory, where the particles have a finite width 
and are modelled by Gaussian densities n(x). For any A 
the corresponding widths are taken from exact diagonal- 
ization calculations yielding = 1 + 2Ay"(|(io|), with 
the mean-field potential 
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n{xi)n{x2) 
Xi~ X2 + \d\ ' 



(8) 



and the pertinent equilibrium distance rfp- Fig. [3] shows 
that this approximation works well for A > 30 (dashed 
line). 

One advantage of these approximations is that they 
can be straightforwardly extended to larger particle num- 
bers. Yet the most interesting parameter range where 
quantum and correlation effects are strong simultane- 
ously (around A = 1) are not accessible by perturbation 
approximations and require a full numerical treatment. 



IV. INFLUENCE OF THE SPIN STATISTICS 
ON THE BREATHING MODE 

A. One-dimensional system 

Let us now repeat the above calculations with excita- 
tions (I) and (II), now starting from a symmetric initial 
coordinate wave function ^fp. According to the Bose- 
Fermi mapping in lcp2j it is expected, that both ini- 
tial states, ^'q and ^f^, should lead to the same QBM. 
However, it is interesting to note that both numerical 
approaches to solve the TDSE (see Sec. II B I fail to 



reproduce this behavior, due to the singularity of the 
Coulomb potential. The standard "regularization" pro- 
cedure which introduces a small finite cut-off n in the po- 
tential, A/ a/ {xi — 2:2)^ -I- K^, yields qualitatively wrong 
results for u!r, cf. Fig. [4] The breathing frequency ojr 
exhibits a nonmonotonic dependence on A, in contrast 
to the solution of the TDSE for an anti-symmetric initial 
state. Reducing the cut-off parameter k slightly improves 
the behavior by reducing the amplitude of the oscillation 
and shifting it towards smaller values of A. But even 
for K as small as 10"^ the spurious oscillation persists 
for A < 1, and no convergence to the behavior of the 
anti-symmetric state is observed. 

We underline that this is not a numerical error and 
not due to the solution procedure, but it is a property of 
the regularized Coulomb potential in ID. The reason for 
the observed unexpected behavior is the incorrect (finite) 
value of the relative wave function (p^ at the origin (i.e. 
at zero particle separation) . This red uces the particle 
repulsion and, thus, ujr, cf. inset of Fig. 



D. 



for more than two particles 



Let us now discuss the dependence of tUr on the number 
of particles N. In the case of classical charged particles 
in a harmonic trap it is well known that the breathing 
frequency is universal, i. e. lo-^ does not depend on N. 
This case describes the limit of strong coupling, A — > 00. 
On the other hand, for the quantum case at finite A, this 
question is still open. It turns out that a numerical so- 
lution of the full time-dependent iV-particle Schrodinger 



B. 2D and 3D systems 

The situation changes completely in 2D and 3D. Now 
particles can avoid each other, allowing for finite values 
foiO), whereas faiO) remains equal to zero. Conse- 
quently, we expect a lowering of uJr for symmetric states, 
compared to anti-symmetric ones, for arbitrary A. This is 
fully confirmed by numerical solutions of the TDSlP^l^ cf. 
Fig. [5j The overall behavior is the same as in ID — we 
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FIG. 5: (color online) QBM frequency Ur, for iV = 2 in a 2D 
trap for symmetric and antisymmetric states from solution 
of the TDSE ([sjl with k = 0. For comparison, also the ID 
result and the two fits i(;/**(A), Eq. are shown by the (red) 
dotted lines. 



observe a QBM with uj]^ — 2 and a second mode Wr(A), 
however, its value is different from ID: the frequency 
is reduced (increased) for an (anti-)symmetric state, cf. 
Fig. [s] While the differences vanish in the classical (and 
quantum) limits due to missing (complete) wave function 
overlap, at intermediate couplings, around A = 1, the 
differences between ID and 2D, as well as between anti- 
symmetric and symmetric states are substantial, reaching 
values of about 3% and 5%, respectively. These differ- 
ences render the QBM frequency uJrW a sensitive di- 
agnostics of the spin statistics of the particles. Indeed, 
in case of a symmetric spin wave function (e.g. spin- 
polarized system in a strong magnetic field), the anti- 
symmetric (symmetric) relative coordinate wave function 
ipo refers to fermions (bosons) and, vice versa, in case of 
an anti-symmetric spin wave function. Furthermore, in 
case of a mixture of fermions and bosons, the resonance 
absorption (II) yields information about the fraction of 
the different components. 

While our results have been obtained for systems with 
a spin-independent hamiltonian, it is straightforward to 
perform analogous TDSE calculations with spin effects 
included. Another interesting observation is that uj,.(X) 
in ID and 2D for the wave function is well described 
by the functional form 

ujI^\\) = aexp[-arctan(6A + c)]-|-d, (9) 
with d = (2 - \/34)/(l - 4), 
dc — exp(7r/2 — arctan c), 
a — (\/3 — d) exp(7r/2), 

cf. Fig.ip. 

Finally, we expect that the reported collective behavior 
will be observable also in anisotropic harmonic traps in d 
dimensions with the difference that there will be a total 
of 2d QBMs with frequencies and w*(A) where i — 
l,...d. 



V. SUMMARY AND DISCUSSION 

In this paper we have presented a complete analysis 
of the quantum breathing modes of interacting charged 
particles in a ID and 2D harmonic trap for arbitrary val- 
ues of the coupling strength A. The results are based on 
accurate solutions of the TDSE with two independent ex- 
citation scenarios and differ qualitatively from previous 
results based on hydrodynamic models. We have shown 
that the frequency ratio 7(A) = ujr{X) /loh, while being 
independent of the trap frequency and the confining sys- 
tem, does depend on the trap dimensionality and on fun- 
damental properties of the particles: their spin statistics 
and the form of their pair interaction. Furthermore the 
A-dependence of 7 is a sensitive measure of the strength 
of many-body effects in the particle ensemble. 

The present results have been obtained for two parti- 
cles with Coulomb interaction where exact solutions of 
the TDSE are possible. Nevertheless, they are represen- 
tative for any nonideal quantum system in a harmonic 
trap. An analogous but different A— dependence exists 
for dipole interaction^^ and for short-range potentials. 
Thus, for a given value of A, Ur is indicative of the form 
of the pair interaction. 

While our results are of importance for any nonideal 
quantum system in a harmonic trap, the y a re of partic- 
ular relevance for electrons^ and excitona^ quantum 
dots as well as for Bose and Fermi systems and their mix- 
tures in traps or optical lattices. There the resonance ab- 
sorption by the QBM, in particular, the relative spectral 
weight of the two modes, may serve as a valuable exper- 
imental diagnostics. Thereby one can take advantage of 
the comparatively easy excitation of the breathing mode 
using the excitation scenarios I or II. 

At the same time the presented precise values of 7(A) 
provide a strong benchmark for nonequilibrium theoreti- 
cal models and computational many-body methods for fi- 
nite systems, including hydrodynamics, TDHF and time- 
dependent densit y fun ctional theory and nonequilibrium 
Greens functions^DUl Indeed, the normal modes of 
trapped systems are an important criterion for the qual- 
ity of nonequilibrium simulations which treat the inter- 
action approximately. This is in complete analogy to ap- 
proximate simulations of atoms which should reproduce 
the excitation energy spectrum as accurately as possible. 
Recently it was shown that the sloshing (Kohn) mode of 
a harmonically trapped system is preserved exactly by 
any conserving many-body approximatiorPSl. Similarly, 
the breathing mode may serve as such a constraint. For 
example, as was shown in Fig. |2] the time-dependent 
Hartree-Fock approximation behaves very well for an 
anti-symmetric initial state (fermions). For symmetric 
states (bosons) Hartree-Fock simulations of strongly cor- 
related charged systems in a harmonic trap have recently 
been developecpS but the results for the breathing mode 
appear to be significantly less accurate than for fermions. 

Finally, a particularly remarkable observation made in 
our paper is the persistence of the center of mass quan- 
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turn breathing mode with frequency lor up to very large 
A values. In strongly coupled systems this mode shows up 
in a "breathing" of each individual particle (see lower left 
part of Fig. [2] and the accompanying videJ^^. Despite 
the fact that for strong coupling, the collective quantum 
degeneracy effects are weak (the interparticle distance 
is much smaller than Ig) each individual particle clearly 
"remembers" its quantum nature. 
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